home *** CD-ROM | disk | FTP | other *** search
/ Language/OS - Multiplatform Resource Library / LANGUAGE OS.iso / dsp / dspgroup / asms.arc / HILBERT.ASM < prev    next >
Encoding:
Assembly Source File  |  1987-04-18  |  6.1 KB  |  387 lines
  1. ;************************************************************
  2. ;             FILTER STRUCTURE
  3. ;
  4. ;
  5. ;
  6. ;             FINITE IMPULSE RESPONSE (FIR)
  7. ;          LINEAR PHASE DIGITAL FILTER DESIGN
  8. ;               REMEZ EXCHANGE ALGORITHM
  9. ;
  10. ;                 HILBERT TRANSFORMER
  11. ;
  12. ;                  FILTER LENGTH =  50
  13. ;
  14. ;             ***** IMPULSE RESPONSE *****
  15. ;           H( 1) =   .13162560E+00 = -H( 50)
  16. ;           H( 2) =   .89106710E-01 = -H( 49)
  17. ;           H( 3) =  -.13654400E+00 = -H( 48)
  18. ;           H( 4) =  -.14179820E-01 = -H( 47)
  19. ;           H( 5) =  -.44728550E-01 = -H( 46)
  20. ;           H( 6) =  -.55178160E-01 = -H( 45)
  21. ;           H( 7) =  -.42822250E-02 = -H( 44)
  22. ;           H( 8) =  -.63906710E-01 = -H( 43)
  23. ;           H( 9) =   .96562650E-02 = -H( 42)
  24. ;           H(10) =  -.56198930E-01 = -H( 41)
  25. ;           H(11) =   .10118810E-01 = -H( 40)
  26. ;           H(12) =  -.38664480E-01 = -H( 39)
  27. ;           H(13) =   .36292080E-02 = -H( 38)
  28. ;           H(14) =  -.13893960E-01 = -H( 37)
  29. ;           H(15) =  -.62213840E-02 = -H( 36)
  30. ;           H(16) =   .17244380E-01 = -H( 35)
  31. ;           H(17) =  -.16327770E-01 = -H( 34)
  32. ;           H(18) =   .54069490E-01 = -H( 33)
  33. ;           H(19) =  -.21842140E-01 = -H( 32)
  34. ;           H(20) =   .97363050E-01 = -H( 31)
  35. ;           H(21) =  -.14074010E-01 = -H( 30)
  36. ;           H(22) =   .15518270E+00 = -H( 29)
  37. ;           H(23) =   .34628240E-01 = -H( 28)
  38. ;           H(24) =   .29468500E+00 = -H( 27)
  39. ;           H(25) =   .54447590E+00 = -H( 26)
  40. ;
  41. ;                       BAND  1       BAND  2       BAND  3       
  42. ; LOWER BAND EDGE      .0000000      .0290000      .4600000
  43. ; UPPER BAND EDGE      .0200000      .4500000      .5000000
  44. ; DESIRED VALUE        .0000000     1.0000000      .0000000
  45. ; WEIGHTING          10.0000000     1.0000000    10.0000000
  46. ; DEVIATION            .0489696      .4896959      .0489696
  47. ;
  48. ; EXTREMAL FREQUENCIES--MAXIMA OF THE ERROR CURVE
  49. ;     .0116667    .0200000    .0290000    .0465000    .0681666
  50. ;     .0898333    .1106666    .1323332    .1531665    .1739999
  51. ;     .1948332    .2156665    .2364998    .2573331    .2781664
  52. ;     .2998331    .3206664    .3414997    .3631664    .3856663
  53. ;     .4073330    .4314996    .4500000    .4600000    .4691667
  54. ;     .5000000
  55. ;
  56. ;***********************************************************
  57. ;
  58. ;
  59. XN:    EQU 0
  60. XNM1:    EQU 1
  61. XNM2:    EQU 2
  62. XNM3:    EQU 3
  63. XNM4:    EQU 4
  64. XNM5:    EQU 5
  65. XNM6:    EQU 6
  66. XNM7:    EQU 7
  67. XNM8:    EQU 8
  68. XNM9:    EQU 9
  69. XNM10:    EQU 10
  70. XNM11:    EQU 11
  71. XNM12:    EQU 12
  72. XNM13:    EQU 13
  73. XNM14:    EQU 14
  74. XNM15:    EQU 15
  75. XNM16:    EQU 16
  76. XNM17:    EQU 17
  77. XNM18:    EQU 18
  78. XNM19:    EQU 19
  79. XNM20:    EQU 20
  80. XNM21:    EQU 21
  81. XNM22:    EQU 22
  82. XNM23:    EQU 23
  83. XNM24:    EQU 24
  84. XNM25:    EQU 25
  85. XNM26:    EQU 26
  86. XNM27:    EQU 27
  87. XNM28:    EQU 28
  88. XNM29:    EQU 29
  89. XNM30:    EQU 30
  90. XNM31:    EQU 31
  91. XNM32:    EQU 32
  92. XNM33:    EQU 33
  93. XNM34:    EQU 34
  94. XNM35:    EQU 35
  95. XNM36:    EQU 36
  96. XNM37:    EQU 37
  97. XNM38:    EQU 38
  98. XNM39:    EQU 39
  99. XNM40:    EQU 40
  100. XNM41:    EQU 41
  101. XNM42:    EQU 42
  102. XNM43:    EQU 43
  103. XNM44:    EQU 44
  104. XNM45:    EQU 45
  105. XNM46:    EQU 46
  106. XNM47:    EQU 47
  107. XNM48:    EQU 48
  108. XNM49:    EQU 49
  109. ;
  110. H0:    EQU 50
  111. H1:    EQU 51
  112. H2:    EQU 52
  113. H3:    EQU 53
  114. H4:    EQU 54
  115. H5:    EQU 55
  116. H6:    EQU 56
  117. H7:    EQU 57
  118. H8:    EQU 58
  119. H9:    EQU 59
  120. H10:    EQU 60
  121. H11:    EQU 61
  122. H12:    EQU 62
  123. H13:    EQU 63
  124. H14:    EQU 64
  125. H15:    EQU 65
  126. H16:    EQU 66
  127. H17:    EQU 67
  128. H18:    EQU 68
  129. H19:    EQU 69
  130. H20:    EQU 70
  131. H21:    EQU 71
  132. H22:    EQU 72
  133. H23:    EQU 73
  134. H24:    EQU 74
  135. ;
  136. MODE:    EQU 75
  137. CLOCK:    EQU 76
  138. MASK1:    EQU 77
  139. MASK2:    EQU 78
  140. TEMPHI:    EQU 79
  141. TEMPLO:    EQU 80
  142. YN:    EQU 89
  143. ONE:    EQU 90
  144. ;
  145.     ORG 0
  146. ;
  147.     B START
  148. ;
  149. ; COEFFICIENTS ARE INITIALLY ;
  150. ; STORED IN PROGRAM MEMORY   ;
  151. ;
  152. ; DUE TO THE SYMMETRY OF THE IMPULSE RESPONSE ;
  153. ; ONLY HALF OF THE SAMPLES OF THE IMPULSE     ;
  154. ; RESPONSE ARE STORED.    THIS MEANS THAT       ;
  155. ; h(N-1-n) = h(n).                  ;
  156. ; 
  157. COEF:
  158.     dw    -4313;            H( 1) =   .13162560E+00 = -H( 50)
  159.     dw    -2919;            H( 2) =   .89106710E-01 = -H( 49)
  160.     dw    4474
  161.     dw    465
  162.     dw    1466
  163.     dw    1808
  164.     dw    140
  165.     dw    2094
  166.     dw    -316
  167.     dw    1841
  168.     dw    -332
  169.     dw    1267
  170.     dw    -119
  171.     dw    455
  172.     dw    204
  173.     dw    -565
  174.     dw    535
  175.     dw    -1771
  176.     dw    716
  177.     dw    -3190
  178.     dw    461
  179.     dw    -5085
  180.     dw    -1135
  181.     dw    -9656
  182.     dw    -17841
  183.     dw    000AH
  184.     dw     0200H
  185.     dw    0800H
  186. ENDC:    dw    08000H
  187. ;
  188.  
  189. START:    LDPK 0
  190. ;
  191.     LACK 1
  192.     SACL ONE    ; CONTENT OF ONE IS 1 ;
  193. ;
  194.     LARK AR0,MASK2    ; THIS SECTION OF CODE LOADS  ;
  195.     LARK AR1,1DH    ; THE FILTER COEFFICIENTS AND ;
  196.     LACK ENDC    ; OTHER VALUES FROM PROGRAM   ;
  197. LOAD:    LARP AR0    ; MEMORY TO DATA MEMORY       ;
  198.     TBLR *-,AR1
  199.     SUB ONE
  200.     BANZ LOAD
  201. ;
  202.     NOP
  203.     NOP
  204. ;
  205. WAIT:    BIOZ NXTPT    ; BIO PIN GOES LOW WHEN A ;
  206.     B WAIT        ; NEW SAMPLE IS AVAILABLE ;
  207. ;
  208. NXTPT:    IN XN,PA3    ; BRING IN THE NEW SAMPLE XN ;
  209.     LAC XN,4
  210.     SUB ONE,15
  211.     SACL XN,0
  212. ;
  213.     ZAC
  214. ;
  215.     LT XNM49    ; DUE TO SYMMETRY h(0) = -h(79) ;
  216.     MPY H0        ; x(n-79) * -h(79) ;
  217. ;
  218.     LTD XNM48
  219.     MPY H1        ; h(1) = -h(78) ;
  220. ;
  221.     LTD XNM47
  222.     MPY H2
  223. ;
  224.     LTD XNM46
  225.     MPY H3
  226. ;
  227.     LTD XNM45
  228.     MPY H4
  229. ;
  230.     LTD XNM44
  231.     MPY H5
  232. ;
  233.     LTD XNM43
  234.     MPY H6
  235. ;
  236.     LTD XNM42
  237.     MPY H7
  238. ;
  239.     LTD XNM41
  240.     MPY H8
  241. ;
  242.     LTD XNM40
  243.     MPY H9
  244. ;
  245.     LTD XNM39
  246.     MPY H10
  247. ;
  248.     LTD XNM38
  249.     MPY H11
  250. ;
  251.     LTD XNM37
  252.     MPY H12
  253. ;
  254.     LTD XNM36
  255.     MPY H13
  256. ;
  257.     LTD XNM35
  258.     MPY H14
  259. ;
  260.     LTD XNM34
  261.     MPY H15
  262. ;
  263.     LTD XNM33
  264.     MPY H16
  265. ;
  266.     LTD XNM32
  267.     MPY H17
  268. ;
  269.     LTD XNM31
  270.     MPY H18
  271. ;
  272.     LTD XNM30
  273.     MPY H19
  274. ;
  275.     LTD XNM29
  276.     MPY H20
  277. ;
  278.     LTD XNM28
  279.     MPY H21
  280. ;
  281.     LTD XNM27
  282.     MPY H22
  283. ;
  284.     LTD XNM26
  285.     MPY H23
  286. ;
  287.     LTD XNM25
  288.     MPY H24
  289. ;
  290.     APAC
  291.     SACH TEMPHI,0
  292.     SACL TEMPLO,0
  293.     ZAC
  294. ;
  295.     LT XNM24
  296.     MPY H24
  297. ;
  298.     LTD XNM23
  299.     MPY H23
  300. ;
  301.     LTD XNM22
  302.     MPY H22
  303. ;
  304.     LTD XNM21
  305.     MPY H21
  306. ;
  307.     LTD XNM20
  308.     MPY H20
  309. ;
  310.     LTD XNM19
  311.     MPY H19
  312. ;
  313.     LTD XNM18
  314.     MPY H18
  315. ;
  316.     LTD XNM17
  317.     MPY H17
  318. ;
  319.     LTD XNM16
  320.     MPY H16
  321. ;
  322.     LTD XNM15
  323.     MPY H15
  324. ;
  325.     LTD XNM14
  326.     MPY H14
  327. ;
  328.     LTD XNM13
  329.     MPY H13
  330. ;
  331.     LTD XNM12
  332.     MPY H12
  333. ;
  334.     LTD XNM11
  335.     MPY H11
  336. ;
  337.     LTD XNM10
  338.     MPY H10
  339. ;
  340.     LTD XNM9
  341.     MPY H9
  342. ;
  343.     LTD XNM8
  344.     MPY H8
  345. ;
  346.     LTD XNM7
  347.     MPY H7
  348. ;
  349.     LTD XNM6
  350.     MPY H6
  351. ;
  352.     LTD XNM5
  353.     MPY H5
  354. ;
  355.     LTD XNM4
  356.     MPY H4
  357. ;
  358.     LTD XNM3
  359.     MPY H3
  360. ;
  361.     LTD XNM2
  362.     MPY H2
  363. ;
  364.     LTD XNM1
  365.     MPY H1
  366. ;
  367.     LTD XN
  368.     MPY H0
  369. ;
  370. ;
  371.     APAC
  372. ;
  373.     SUBS TEMPLO
  374.     SUBH TEMPHI
  375.     
  376.     
  377.     SACH YN,1
  378.     LAC YN,12
  379.     ADDH MASK1
  380.     SACH YN,0
  381. ;
  382.     OUT YN,PA4    ; OUTPUT THE FILTER RESPONSE y(n) ;
  383. ;
  384.     B WAIT        ; GO GET THE NEXT POINT ;
  385. ;
  386.     END
  387.